New CMB constraints on the cosmic matter budget: trouble for nucleosynthesis? 
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We compute the joint constraints on ten cosmological parameters from the latest CMB measure- 
ments. The lack of a significant second acoustic peak in the latest Boomerang and Maxima data 
favors models with more baryons than Big Bang nucleosynthesis predicts, almost independently 
of what prior information is included. The simplest flat inflation models with purely scalar scale- 
invariant fluctuations prefer a baryon density 0.022 < h?Q,i, < 0.040 and a total nonbaryonic (hot + 
cold) dark matter density 0.14 < h^Qdm < 0.32 at 95% confidence, and allow reionization no earlier 
than 2 ~ 30. 



One of the main challenges in modern cosmology is to 
refine and test the standard model of structure formation 
by precision measurements of its free parameters. The 
cosmic matter budget involves at least the four parame- 
ters rib, f^cdm, f^i/ and r^A, which give the percentages of 
critical density corresponding to baryons, cold dark mat- 
ter, massive neutrinos and vacuum energy. A "budget 



deficit" 



1 — — Hp 



n^, — i}\ manifests itself 



as spatial curvature. The description of the initial seed 
fluctuations predicted by inflation requires at least four 
parameters, the amplitudes As & At and slopes Ug & nt 
of scalar and tensor fluctuations, respectively. Finally, 
the optical depth parameter r quantifies when the first 
stars or quasars reionized the Universe and the Hubble 
parameter h gives its current expansion rate. 
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The 87 band power measurements used, 
shows the simple inflationary model with 
■ = 0, = 0.43, h^n^dm = 0.20, 
1, h — 0.63. Note that although we 



include the calibration uncertainties in our analysis, they are 
not reflected by the plotted error bars. 



During the past year or so, a number of papers 
have used the measured cosmic microwave background 
(CMB) fluctuations to constrain subsets of these param- 
eters. CMB data has improved dramatically since fluctu- 
ations were first detected . The measurement of a first 
acoustic peak at the degree scale , suggesting that the 
Universe is flat (ilk = 0), has now been beautifully con- 
firmed and improved by using the ground-breaking high 
fidelity maps of the Boomerang and Maxima jlj] 
experiments. As can be seen in Figure |^, perhaps the 
most important new information from Boomerang and 
Maxima is their accurate measurements of the angular 
power spectrum Ci on even smaller scales, out to mul- 
tipole £ ^ 600 — 800. The striking lack of a significant 
second acoustic peak places strong constraints on the cos- 
mological parameters , making a new full-fledged 
analysis of all the CMB data very timely. 

In this Letter, we jointly constrain the following 10 cos- 
mological parameters: r, ilk, ^a, ns, rit, As, the tensor- 
to scalar ratio r = At /As, and the physical matter den- 
sities = h^i^h, Wcdm = h'^^cdm and LOi, = h'^fl^. The 
identity h = \/ (wcdm + Wb + Wi/)/(l — — ^a) fixes 
the Hubble parameter. We use the 10-dimensional grid 
method described in In essence, this utilizes a tech- 
nique for accelerating the the CMBfast package |Q by a 
factor around 10'^ to compute theoretical power spectra 
on a grid in the 10-dimensional parameter space, fitting 
these models to the data and then using cubic interpo- 
lation of the resulting 10-dimensional likelihood function 
to marginalize it down to constraints on individual or 
pairs of parameters. We use the 87 data points shown 
in Figure |l|, combining the 65 tabulated in |^] with the 
12 new Boomerang points ijlll and the 10 new Maxima 
points |]l2|. 

Our 95% confidence limits on the best constrained pa- 
rameters are summarized in Table 1. Figure ^ shows 
that CMB alone suggests that the Universe is either flat 
(near the diagonal line + ^A = 1, where f2m = 
ilb + i^cdm + ^ly) or closed (upper right). These con- 
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straints come largely from the location of the first peak, 
which is well-known to move to the right if the curvature 
is increased |l^. Very closed models work only be- 
cause the first acoustic peak can also be moved to the 
right by increasing the tilt Us or decreasing the matter 
density and bringing the large-scale COBE signal back up 
with tensor fluctuations (gravity waves) [|l6|,^ . Galaxy 
clustering constraints disfavor such strong blue-tilting, 
and Figures ^ and || show that closing this loophole by 
barring gravity waves (r — 0) favors curvature near zero 
and Us near unity. This is a striking success for the old- 
est and simplest inflation models, which make the three 
predictions r « 0, Jlk ~ and rig « 1 |l^,|l^. Another 
important success for inflation is that the first peak is 
so narrow — if the data had revealed the type of broad 
peak expected in many topological defect scenarios, none 
of the models in our grid would have provided an accept- 
able fit. Because of these tantalizing hints that "back to 
basics" inflation is correct. Table 1 and Figure || include 
results assuming this inflation prior r = U,]^ — Q,ns = 1. 

Table 1 - Maximum-likelihood values and 95% confidence limits. 
The "inflation prior" for each parameter is indicated in boxes in 
Figure ^. Vldm = f^cdm -I- A dash indicates that no limit was 
found, with the likelihood still above at the edge of our grid. 
Extrapolation would suggest a limit ris ^ 1.75. 



Quantity 


10 free parameters 
Min Best Max 


Inflation prior 
Min Best Max 


T 


0.0 0.0 0.33 
.017 .05 
0.02 0.08 
- 0.2 0.80 


0.0 0.0 0.28 
.022 .03 .040 
0.14 0.20 0.32 
-0.16 .43 0.65 




- -0.6 0.13 
0.8 1.5 


-0.13 0.10 
0.84 1.0 1.17 



The constraints in Table 1 are seen to be much more 
interesting than those before Boomerang and Maxima 
1^, thanks to new information on the scale of the sec- 
ond peak and beyond. Cold dark matter and neu- 
trinos have indistinguishable effects on the CMB ex- 
cept for very light neutrinos (small w^), and the cur- 
rent data still lacks the precision to detect this subtle 
difference. The predicted height ratio of the first two 
peaks therefore depends essentially on only three pa- 
rameters [p|, p^p3| : Us, tOh and uJdm, where the total 
dark matter density ujdrn = ^^cdm + w^. Let us focus 
on the constraints on these parameters. Increasing Wb 
tends to boost the odd-numbered peaks (1, 3, etc.) at 
the expense of even ones (2, 4, etc.) whereas in- 

creasing codm suppresses all peaks (see the CMB movies 
at www.hep.upenn.edu/ '^max or www.ias.edu/ ^whu). 
The low second peak can therefore be fit by either de- 
creasing the tilt Us or by increasing the baryon density 
I^Jl^ compared to the usually assumed values « 1, 
Wb ~ 0.02. As illustrated in Figure ^ this conclusion is 



essentially independent of what priors are assumed. How- 
ever, reducing rig below 0.9 is seen to make things worse 
again, as the first peak becomes too low relative to the 
COBE-normalization. 




FIG. 2. The regions in the (Jim, J^A)-plane that are ruled 
out at 95% are shown using (starting from the outside) no 
priors, the prior that 0.5 < h < 0.8 (95%), and the additional 
constraint r = The SN la constraints are from White [20] 
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FIG. 3. Marginalized likelihoods assuming that 

0.5 < h < 0.8 (95%) and the inflationary priors specified 
in the boxes. 2 — a limits are roughly where the curves drop 
beneath the dashed lines. 

In short, there are two very simple ways of explain- 
ing the lack of a prominent second acoustic peak: more 
baryons or a red-tilted spectrum |P,p^. However, as we 
will now discuss, both of these solutions have problems 
of their own. 
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Figure |5| shows that when more baryons are added, 
more dark matter is needed to keep the first peak height 
constant. When the tilt Ug is fixed by the inflation prior, 
the constraints on the remaining two parameters tOdm 
and (jJh are seen to become quite tight. Intriguingly, 
the preferred baryon fraction is of the same order as 
preferred by Big Bang nucleosynthesis, but nonetheless 
higher than the tight nucleosynthesis error bars [|^j2| 
tJb = 0.019 ± 0.0024 allow. Even if the nucleosynthesis 
error bars have somehow been underestimated so that 
Wb ^ 0.023 as required by the CMB data plus simple 
inflation is allowed, this solution may conflict with other 
astrophysical constraints. For instance, X-ray observa- 
tions of clusters of galaxies can be used to determine the 
ratio of baryons to dark matter [^,|ol, and uJh = 0.03 
can only be reconciled with these observations by having 
Qm ^ 0.7 which would conflict with the supernova la 
results and other estimates of the dark matter density 

ill- 

On the other hand, the tilt solution is no panacea ei- 
ther. In a class of popular inflationary models known 
as power law inflation, the amplitude of the tensor com- 
ponent is approximately related to the tilt of the scalar 
spectrum, r ~ 7(1 — rig) If we choose to fit the 

data by lowering the tilt to Ug — 0.9, this would raise the 
COBE-normalization by 70%. Models that match the 
COBE normalization therefore make the first peak too 
low by a factor of 1.7 in power, which is ruled out by the 
data. In other words, imposing r ~ 7(1 — n.,) (which we 
have not done in our analysis) would exclude as low 
as 0.02. Thus the simple tilt solution does not work for 
all inflation models. 

Could the apparent problem be a mere statistical 
fluke? It would certainly be premature to claim a rock- 
solid discrepancy between CMB and nucleosynthesis plus 
power law inflation. The x^-value for the best flt infla- 
tion model with Wb = 0.02 is still statistically acceptable 
« 81 for 87 degrees of freedom reduced by about 
5 effective parameters). However, serious discrepancies 
in peak heights tend to get statistically diluted by the 
swarm of points with large error bars at lower £ that agree 
with most anything reasonable (indeed, drops down 
to 71 for Wb — 0.03 and as low as 68 without any priors), 
and the relative likelihood rises sharply with uJh regard- 
less of what priors are imposed. To assess the sensitivity 
of the results to the choice of data, we therefore repeated 
our entire analysis for the following cases: (a) using all 
the data except Maxima and (b) using only COBE and 
the new Boomerang data. Omitting Maxima removed 
the "CMB only" and "CMB+/i" exclusion regions that 
are seen to protrude in from the left in Figure ^. This 
is because the Maxima points place an upper limit on 
the height of (the left part of) the third peak, effectively 
giving an upper limit on the baryon density. Dropping 
Maxima also loosened the upper limit on ujdm somewhat 
and marginally weakened the bounds on f2k and Ha ■ The 



other constraints were essentially unaffected. Most im- 
portantly, the lower bound on Wb seen in Figure ^ re- 
mained unchanged, since it comes from the low ratio of 
the 2nd to 1st peak heights |23) . 




Scalar spectral index 



FIG. 4. The regions in the (ris, (jjb)-plane that are ruled 
out at 95% are shown using (starting from the outside) no 
priors, the prior 0.5 < h < 0.8 (95%), the additional con- 
straint r = and the additional constraints r = (dashed 
line). The horizontal band shows the nucleosynthesis con- 
straints t^b = 0.019 ± 0.0024. 

Although our inclusion of the 10% uncertainly in the 
Boomerang's calibration (20% in power) was not very 
important in our full analysis, as the fitting procedure 
de facto calibrated Boomerang off of other experiments, 
this substantially degraded the results of our COBE -|- 
Boomerang analysis. We therefore repeated it three more 
times, without the calibration error but multiplying the 
Boomerang points by 0.9, 1.0 and 1.1, respectively. The 
results were quite similar to those using all the data, as 
expected from the experimental concordance seen in Fig- 
ure ^ However, most constraints got slightly tighter, 
consistent with the above-mentioned dilution hypoth- 
esis. Rather than go away, the baryon problem became 
exacerbated: the 95% inflationary lower limit on lu^, was 
tightened from 0.024 to 0.027 with x^ = 12 (with a total 
of 20 Boomerang -I- COBE points and 4 free parameters). 
In contrast, the tilt solution gave = 22, and higher still 
when the Boomerang normalization was raised or lowered 
by 10%. Most strikingly, in the COBE+Boomerang ver- 
sion of Figure ^, tJb is not permitted to be low enough 
to agree with nucleosynthesis for any value of the tilt rig , 
so the tilt solution may have worked using all the data 
merely because of the above-mentioned dilution effect. 

Can the baryon problem be explained by inaccuracies 
in our numerical method? The correlations between the 
Boomerang points (which we could not include since the 
have not yet been made public) are reportedly very small 
[prf. Although a range of approximations are involved 
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as detailed in g, for instance in the likelihood calcula- 
tion, it appears unlikely that such inaccuracies are large 
enough to have a major impact on the lower bound on 
ujh- Perhaps the best indication of this is that a number 
of independent analyses pl|-p^ have been made available 
since this paper was originally submitted, using a wide 
range of computational techniques, and they all favor 
baryon fractions in excess of the current nucleosynthesis 
prediction. More baryons also solve some older problems 

i5i. 
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FIG. 5. 
(tJdm, t^b)-plane. 
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Dark matter density h 
Similar to Figure 4, but for for the 



A number of other ways out have been proposed |H,|T3[ , 
ranging from mechanisms for delaying standard recom- 
bination to time-variation of physical constants and se- 
vere mis-estimates of the Boomerang beam width. How- 
ever, these explanations are all of a highly speculative 
nature. An excellent way to clear up this mystery will 
be to search for a third acoustic peak, which is boosted 
by more baryons but suppressed by most of the other 
proposed remedies. 

Apart from the matter budget. Table 1 and Figure ^ 
also show that the CMB data provides perhaps the first 
meaningful upper limit on r, the optical depth due to 
reionization (compare |||,^). Since r oc /iSlt-^km ^m"^^^ 
if the redshift of reionization Zjon 3> 1, our lower limit 
ujh > 0.024 combined with our upper limit r < 0.35 
give the constraint Zion ^ Adh'^^^H.U^ , or Zjon ^ 28 for 
h = JIa = 0.7. This is compatible with the range Zjon = 
8 — 20 favored by numerical simulations, but challenges 
more extreme models. 

In conclusion, the new Boomerang results look like a 
triumph for the simplest possible inflationary model but 
for one rather large fly in the ointment: the lack of a 
significant second acoustic peak suggests that we may 
need to abandon either a popular version of inflation, 
the current nucleosynthesis constraints, or some even 
more cherished assumption. In answering one question. 
Boomerang has raised another. Its answer is likely to lie 
in the third peak, and the race to reach it has now begun. 
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